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Abstract 

Using an isospin- and momentum-dependent transport model we study ef- 
fects of the momentum-dependent symmetry potential on heavy-ion collisions 
^ ' induced by neutron-rich nuclei. It is found that symmetry potentials with and 

c ■ 

without the momentum-dependence but corresponding to the same density- 
^ \ dependent symmetry energy Esym{p) lead to significantly different predictions 

on several i?sj/m(/o)-sensitive experimental observables especially for energetic 
nucleons. The momentum- and density-dependence of the symmetry poten- 
tial have to be determined simultaneously in order to extract the Esym{p) 
accurately. The isospin asymmetry of midrapidity nucleons at high trans- 
verse momenta is particularly sensitive to the momentum-dependence of the 
symmetry potential. It is thus very useful for investigating accurately the 
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equation of state of dense neutron-rich matter. 
25.70.-z, 25.70.Pq., 24.10.Lx 
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I. INTRODUCTION 



The equation of state (EOS) of asymmetric nuclear matter at density p and isospin 
asymmetry 5 = {Pn — Pp) / {Pp + Pn) is usually expressed as 

E{p, 5) = E{p, 5 = 0) + E,,M5' + 0{5^), (1) 

where £^sym(p) is the nuclear symmetry energy. The latter is very important for many in- 
teresting astrophysical problems and yet poorly known for dense neutron-rich matter [1-3] . 
Predictions based on various many-body theories diverge widely, especially at high densities. 
In fact, even the sign of the symmetry energy above 3po is still very uncertain [4]. Since 
nuclear reactions induced by high energy radioactive beams can produce transiently dense 
neutron-rich matter, the planned Rare Isotope Accelerator (RIA) and the new accelerator 
facility at GSI provide the first opportunity in terrestrial laboratories to explore experi- 
mentally the EOS of dense neutron-rich matter [5-8]. Crucial to the extraction of critical 
information about the £'sym(p) is the comparison between experimental data and transport 
model calculations. Based on isospin-dcpendent transport model calculations, several ex- 
perimental observables have been identified as promising probes of the Egymip), such as, the 
neutron/proton ratio [9], isoscaling in nuclar multifragmentation [10-13], the neutron-proton 
differential flow [14-16], the neutron-proton correlation function [17] and the isobaric yield 
ratios of light clusters [18]. 

An important input to the transport models are the isovector (symmetry) and isoscalar 
potentials. However, in all transport models for heavy-ion collisions the momentum- 
dependence of the symmetry potential was seldom taken into account. In the recent work 
of ref. [19] effects of the neutron- proton effective mass splitting due to the momentum- 
dependent nuclear potentials were studied using Bombaci's parameterization of single nu- 
cleon potential [4]. Nevertheless, effects of the momentum- dependence of the symmetry 
potential in heavy-ion collisions are still largely unknown. On the other hand, effects of the 
momentum-dependence of the isoscalar potential in nuclear reactions are well-known [20-26]. 
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This is mainly because only very recently was the momentum-dependence of the symmetry 
potential given in a form practically usable in transport model calculations [27]. In this work, 
we study in detail effects of the momentum-dependence of the symmetry potential on heavy- 
ion collisions induced by neutron-rich nuclei. A short report of this study can be found in 
ref. [28] . In the next section, we discuss the momentum-dependence of the symmetry poten- 
tial using as a guide results of Hartree-Fock calculations with Gogny effective interactions. 
The strengths of the symmetry potential with and without the momentum-dependence are 
compared. In section III, we present and discuss effects of the momentum-dependence of 
the symmetry potential on several experimental observables in heavy-ion collisions induced 
by neutron-rich nuclei. A summary is given at the end. 



II. MOMENTUM-DEPENDENCE OF THE SYMMETRY POTENTIAL 

That the momentum-dependence of the single particle potential is different for neutrons 
and protons in asymmetric nuclear matter is well-known and has been a subject of intensive 
research based on various many-body theories using non-local interactions, see e.g., ref. [4] 
for a review. Experimentally, it is known that the strength of nuclear symmetry potential 
decreases hnearly from about 28 ± 6 MeV with increasing kinetic energy with a slope of 
about 0.1-0.2 from optical model analyses of nucleon-nucleus scattering data [29]. However, 
the available data is limited only to nucleon energies less than 100 MeV. Guided by a 
Hartree-Fock calculation using the Gogny effective interaction, the single nucleon potential 
was recently parameterized as [27] 

U{p,S,p,T)^A^^ + Ai^ 
Po Po 



2^-' ^dv, (2) 
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The parameter x can be adjusted to mimic predictions on the density dependence of 
the symmetry energy Egymip) by microscopic and/or phenomenlogical many-body theo- 
ries. In the above r = 1/2 (—1/2) for neutrons (protons) and t ^ t'; a = 4/3; JrifiP) 
is the phase space distribution function at coordinate r and momentum p. The parame- 
ters Au,Ai,B, Cr^T, Cry and A were obtained by fitting the momentum-dependence of the 
U{p,5,p,t) predicted by the Gogny Hartree-Fock and/or the Brueckner-Hartree-Fock cal- 
culations, saturation properties of symmetric nuclear matter and the symmetry energy of 
30 MeV at normal nuclear matter density po — 0.16/ fm^ [27]. More specifically, B=106.35 
MeV and A = is the nucleon Fermi momentum in symmetric nuclear matter. Au and Ai 
depend on the x parameter according to 

2iB 

Au{x) = -95.98 - X , (3) 

cr + I 

IB 

AAx) = -120.57 + X . (4) 

The compressibility of symmetric nuclear matter Kq is set to be 211 MeV. 

The momentum-dependence of the symmetry potential steams from the different inter- 
action strength parameters and C^^t for a nucleon of isospin r interacting, respec- 
tively, with unlike and like nucleons in the background fields. More specifically, we obtained 
Cuniike = —103.4 MeV and Cuke = —11.7 MeV. One characteristic of the momentum depen- 
dence of the symmetry potential is the different effective masses for neutrons and protons 
in isospin asymmetric nuclear matter. In terms of the single particle potential U (p, 5, p, r) 
the nucleon effective mass is defined as 

-1 



rrir 



rrir dU 
p dp 



(5) 

Since the momentum dependent part of the single particle potential in eq. 2 is independent 



of the parameter x, the nucleon effective mass is the same for different symmetry energies. 
Shown in Fig. 1 are the nucleon effective masses as a function of density (upper window) 
and isospin asymmetry (lower window). These results are consistent with the Bruckner 
Hartree-Fock predictions [37]. 
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FIG. 1. Nucleon effective masses as a function of density (upper window) and isospin asymme- 
try (lower window). 

The parameter x was introduced to reflect the largely uncertain behavior of Esym{p) at 
supranormal densities, specifically its potential part E^y^{p). The latter corresponding to 
X = 1 (denoted by MDI(l) which gives the same Esym{p) as the default Gogny interaction) 
and X = (MDI(O)) are shown in Fig. 2. For comparisons we have also included the 
potential part of the more stiffer symmetry energy Efy^^{p) = 30p/po, predicted by several 
models, such as, the Relativistic-Mean-Field (RMF) model. As a reference, the kinetic 
contribution Egy^(p) = /i^/6m„(37r^p/2)^/^ is also shown. The variation of Egymip) resulted 
by changing from the MDI(l) to the MDI(O) parameter set is well within the uncertain 
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range predicted by various many-body theories, see e.g., refs. [32,33,37]. We found that the 
potential contribution to the symmetry energy can be parameterized as 



EPy^{p) = 3.08 + 39.6m - 29.2^^ + 5.68m^ - 0.523u^ 



for the MDI(l) (Gogny) interaction and 



EP (p) = 1.27 + 25.4m - 9.31m^ + 2.17u^ - 0.21m^ 



for the MDI(O) interaction, where u = p/po is the reduced nucleon density. 
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FIG. 2. The density dependence of the potential and kinetic parts of nuclear symmetry energy. 



To compare with the experimentally found momentum-dependence of the symmetry 
potential in nucleon-nucleus scatterings, we illustrate in Fig. 3 the strength of the symmetry 
potential {Un — Up)/25 as a function of nucleon kinetic energy at densities of 0.5po) Po and 2po, 
for X = 1 (left) and x = (right), respectively. The results are obtained by using 5 = 0.2. 
However, since the isovector part of the nucleon potential is approximately propotional to 



S, the quantity ([/„ — Up)/ 25 is approximately independent of 5. At normal density, for 
both X = 1 and x = the strength of symmetry potential decreases from about 33 MeV 
at Ekin = with a slope of about 0.22 for Ekm < 100 MeV in general agreement with 
the nucleon-nucleus scattering data. Moreover, it is seen that the strength approaches zero 
above about 200 MeV. At other densities, however, the predictions are dramatically different 
for X = 1 and x = especially at high densities. For p = 2po, for instance, the strength 
of symmetry potential is negative and decreases quickly with energy indicating a stronger 
attractive/repulsive symmetry potentials for neutrons/protons with x = 1. With x = 0, 
however, at the same density p = 2po the strength of the symmetry potential is now positive 
indicating the repulsive/attractive nature of the symmetry potential for neutrons/protons. 
It is thus characteristically different from the x = 1 case. It is also seen that the strength 
of the symmetry potential decreases relatively slowly with x = 0. These different energy 
dependences of the symmetry potential are correlated to the negative and positive potential 
contribution to the symmetry energy at high densities for x = 1 and a; = 0, respectively. 
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FIG. 3. Strength of the symmetry potential as a function of kinetic energy at three densities 
for X = 1 (left) and x = (right). 

Since the C-terms in the single particle potential depend inseparably on the density, 
momentum and isospin, to investigate effects of the momentum-dependence of the symmetry 
potential we shall also compare results obtained using eq. 2 with those obtained using a 
single nucleon potential Unomsip, ^,P, t) = Uo{p,p) + Usymip, ^, that has almost the same 
momentum-dependent isoscalar part Uo{p,p) as that embedded in eq. 2 and a momentum- 
independent symmetry potential Usym{p,S,T) that gives the same Egymip) as eq. 2. The 
momentum-independent Usym{p,S,T) can be obtained from 

Usym{P: t) = dWsym/dpr: (8) 

where Wsym is the isospin-dependent part of the potential energy density 

Wsym^E^^^M-P-S'. (9) 

For the symmetry energy Ef^^{p), for instance, one obtains 

Km^iP^ (5, r) = 4T5(30ii - I2.1u^/'') + 4.231^2/35'. (10) 

Using the parameterizations for EPy^{p) in eqs. 6 and 7, we obtain 

U^m\p, ^) = 4t5(3.08 + 39.6ii - 2%.2u^ + bMu^ 

-0.52m^) - (5^(3.08 + 29.2^2 - WAii^ + l.blu^) (11) 

for the MDI(l) parameter set and 

U^m^P, 5> ^) = 4r5(1.27 + 25.4m - 9.31m' + 2.17u^ 

-0.21m^) - 5^(1.27 + 9.31m' - 4.33^^ + 0.63m^) (12) 

for the MDI(O) parameter set. The momentum-dependent isoscalar potential Uq{p,p) will 
be given in the following shortly. 



9 



The strengths of the symmetry potentials with and without the momentum-dependence 
but corresponding to the same Egymip) are compared in Fig. 4 by examining the difference be- 
tween neutron and proton potentials (f/„ — f/p)/25. Again, this quantity is approximately in- 
dependent of the parameter 5. The symmetry potential without the momentum- dependence 
is higher in magnitude and has generally steeper slopes than the momentum-dependent one 
for p/po ^ 2.3 with the MDI(l) parameter set and at all densities for the MDI(O) parameter 
set. Moreover, the strength of momentum-dependent symmetry potentials decreases with 
the increasing momentum. Thus the difference between the symmetry potentials with and 
without the momentum-dependence is larger for nucleons with higher momenta. 




FIG. 4. .Strengths of the symmetry potentials with and without (solid) the momen- 
tum-dependence as a function of density as measured by the difference between neutron and proton 
potentials. 

To identify reliably effects of the momentum-dependence of the symmetry potential with- 
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out much interference from density effects the Unomsip, SyP, t) should lead to about the same 
reaction dynamics and evolution of density profiles as the single particle potential in eq. 2. 
Both of them are mainly determined by the isoscalar potential for which we select the original 
momentum-dependent Yukawa interaction (MDYI) [22] 



1.24 



Uoip,p) = -110.44m + 140.9m 

po J ^l + (p-p')V(1.58pO^)2' ^ ""^ 

where p% is the Fermi momentum. The compressibility Kq for this interaction is 215 MeV. 
This potential is compared with the variational calculation using the UVIA + UVII inter- 
actions in Fig. 5. Except at high densities but low momenta the MDYI interaction gives 
about the same isoscalar potential as the variational approach as first found in ref. [22]. 
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FIG. 5. Strengths of the isoscalar potential as a function of density at the four wave numbers 
for the MDI(l) and MDI(O) parameter sets (lower window) as well as the MDYI (upper window) 
interaction in comparison with the variational many-body calculations using the UVIA two-body 
and the UVII three-body potentials, (upper window) 

In the lower window of Fig. 5 we compare the isoscalar part (C/„ + Up)/ 2 of the single 
nucleon potentials corresponding to the MDI(l) and MDI(O) parameter sets. It is seen that 
they have roughly the same values at all wave numbers k examined, except at high densities 
where the MDI(O) parameter set leads to slightly stiffer isoscalar potentials. This feature 
ensures that the differences we see in reaction calculations using the MDI(l) and MDI(O) 
parameter sets are truly due to the different symmetry potentials. As a reference we have also 
included Wiringa's calculations of the isoscalar potentials within the variational many-body 
approach using the UVIA two-body potential and the UVII three-body potential. The single 
nucleon potential was obtained by using a realistic Hamiltonian that fits nucleon-nucleon 
scattering data, few-body nuclear binding energies and nuclear matter saturation properties 
[30]. It was shown that combinations of different two-body and three-body forces lead to 
slightly different single nucleon potential especially at high densities and/or high momenta. 
It is seen that for wave numbers less than about 2fm~^ the MDI(l) and MDI(O) results 
are very close to that of the variational calculations. At higher wave numbers, however, 
the variational results are higher. For our studies of isospin effects in this work isoscalar 
potentials are less important as long as the reaction dynamics is about right. We thus 
consider the isoscalar parts of both the MDI(l) and MDI(O) interactions as very reasonable 
for the purpose of this work. 

Overall, the MDYI interaction leads to about the same isoscalar potential as the MDI(l) 
and MDI(O) interactions for k less than about 2fm~^. We shall also verify numerically that 
the potential Unoms{p, ^, P, t) constructed this way indeed leads to about the same reaction 
dynamics and the evolution of density profiles as the potential in eq. 2. For comparisons and 
identify effects due to the momentum-dependence and stiffness of the isoscalar potential, we 
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use occasionally also the momentum-independent Bertsch-Kruse-DasGupta (BKD) isoscalar 
potentials. The soft BKD is given by [31] 



(14) 



and stiff BKD is given by 



(15) 



The compressibihties of these two potentials are Kq — 201 and 377 MeV, respectively. 



In this section we examine effects of the momentum-dependence of the symmetry po- 
tential on the reaction dynamics and several experimental observables. Several observables 
are known to be sensitive only to the symmetry potential but not to the isoscalar poten- 
tial. They are mainly neutron-proton differential or relative quantities [9,14] where effects 
of the isoscalar potential with or without the momentum-dependence are largely canceled 
out. Moreover, these observables are also not sensitive to the in-medium nucleon-nucleon 
cross sections [9,17,18] since the latter act approximately identically on both neutrons and 
protons. Assuming the symmetry potential is momentum- independent, these observables 
were previously proposed as promising probes of the density- dependence of the symmetry 
energy. Our emphasis in this work will be on searching for experimental observables that can 
be used to determine simultaneously both the density- and momentum-dependence of the 
symmetry potential, and thus to determine uniquely the density-dependence of the symme- 
try energy. These studies are carried out within an isospin-dependent transport model [9], 
for the latest review of this model we refer the reader to ref . [7] . In this work we concentrate 
on ^^^Sn-l-^^'^Sn reactions at the highest RIA beam energy of 400 MeV/nucleon. 



III. RESULTS AND DISCUSSIONS 
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A. Effects on reaction dynamics 

Shown in Fig. 6 are the evolutions of central baryon densities in ^^^Sn+^^^Sn reactions at 
a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm with various potentials. 
In the upper window, we first compare results obtained by using the soft BKD (solid thin 
line) and the MDYI (dot-dashed line) isoscalar potentials but the same symmetry potential 
^sym^ip)- It is seen that the soft BKD leads to a higher compression by about 10% compared 
to the MDYI interaction in agreement with previous findings in many studies. 




FIG. 6. Evolution of central baryon density with various potentials in the Sn+ Sn reaction 
at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 

Next, we compare results obtained with the MDI(O) (long dashed line) and the MDI(l) 
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(dashed line) potentials. As we have discussed earher they have about the same isoscalar 
potential while very different symmetry potentials as shown in Fig. 4. It is interesting to 
notice that starting at about 10 fm/c when the central density is above about 2.4po, the 
central density for the MDI(l) starts increasing significantly with respect to the MDI(O) case. 
This is what we have expected based on eq. 1 and the density dependence of the symmetry 
energy. The relative values of the central densities reflect closely the relative stiffnesses of the 
underlying equations of state. The above observation is due to a softening of the underlying 
EOS in the MDI(l) case. As shown in Fig. 2, the potential part of the symmetry energy 
increases continuously for the MDI(O) while decreases with increasing density above po for the 
MDI(l) case. In fact, starting at about 2.4po, the symmetry energy becomes negative with 
the MDI(l) parameter set, leading thus to the more compressed central region. Overall, 
however, the differences in central densities using different potentials are less than about 
10%. 

In the lower window of Fig. 6, we study effects of the momentum-dependence of the 
symmetry potential on the central density. This is done by comparing results obtained 
by using the MDI(O) with MDYI+U^^^^°\ and MDI(l) with MDYI+C/j^^^^^), respectively. 
The two pairs of potentials correspond to the same symmetry energy of MDI(O) and MDI(l), 
respectively. It is seen that without the momentum-dependence of the symmetry potential, 
i.e., MDYI+C/j^^^(^) and MDYI+C/j^^^(°) cases, the central density is slightly lower. This is 
because the symmetry potential is effectively stronger without the momentum-dependence 
as shown in Fig. 4. For neutron-rich matter this effectively stiffens the EOS, leading therefore 
to a lower compression. However, the general dynamics of the reaction is merely affected by 
the momentum-dependence of the symmetry potential. 

B. Effects on nucleon emissions 

The generally repulsive/ attractive symmetry potential for neutrons/protons at subnor- 
mal densities is to cause more neutrons/protons to be free/bound. We identify free nucleons 
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as those with local densities less than po/8. At midrapidity, the isospin asymmetry of free 
nucleons are dominated by the high-density behavior of nuclear symmetry energy. As the 
symmetry energy becomes more stiffer from the MDI(l) through MDI(O) to the Ef^^[p) as 
shown in Fig. 2, gradually more/less neutrons/protons are expected to be emitted at midra- 
pidity. The different symmetry energies will thus lead to quite different isospin asymmetries 
of free nucleons at mid-rapidity. Shown in the upper window of Fig. 7 are the nucleon 
rapidity distributions with the four different potentials. Since the three potentials MDI(l), 
MDI(O) and MDYI + U^'^^{p) have about the same isoscalar part while their isovector 
parts are quite different, the observed difference among results using them clearly reveals 
the role of the symmetry potential in nucleon emissions. We also included the calculation 
with the SBKD -\- U^^^{p) interaction as a reference. First of all, it is seen that midra- 
pidity nucleons are most sensitive to the variation of symmetry energy as one expects since 
these nucleons are the ones having gone through the high density phase. Moreover, neutron 
rapidity distributions are more affected by the symmetry energy than protons. The dN/dy 
values around midrapidity are higher for neutrons with the stiffer symmetry energies as one 
expects. 
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FIG. 7. Rapidity distributions of neutrons and protons with various potentials in the 
^^^^n -|-124 reaction at a beam energy of 400 MeV/nucleon and an impact parameter of 1 
fm. 

Included in the upper window Fig. 7 are also results obtained with the soft BKD isoscalar 
potential and the U^^^{p) symmetry potential. Compared to the calculations with the 
MDYI isoscalar potential but the same symmetry potential U^jl^^{p)i the soft BKD leads 
to reduced emissions of both neutrons and protons. In particular, the resulting neutron 
rapidity distribution is very close to the MDI(O) case. This illustrates one serious drawback 
of inclusive observabes for studying the isospin-dependence of the nuclear EOS, such as the 
nucleon rapidity distributions studied here, that they are sensitive to both the isoscalar and 
isovector potentials. This is an important reason for us to concentrate on observables that 
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are only sensitive to the symmetry potential in the following sections. 

What are the effects of the momentum-dependence of the symmetry potential on nu- 
cleon rapidity distributions? To answer this question we compare in the lower window 
of Fig. 7 nucleon rapidity distributions obtained with (MDI(O) and MDI(l)) and without 
(MDYI + U^^^^^^ and MDYI + U^^^^^^^) the momentum-dependence in the symmetry 
potentials. It is seen that the momentum-dependence of the symmetry potential has very 
little effect on protons. This is mainly because the Coulomb potential works against and 
dominates over the symmetry potential on protons. However, the momentum-dependence 
of the symmetry potential affects significantly the neutron rapidity distributions. Overall, 
without the momentum-dependence of the symmetry potential (solid and dotted lines) the 
dN/Dy values for neutrons increases by 10 to 15% with respect to the results with the 
momentum-dependent symmetry potentials (long dashed and dashed lines). This increased 
emissions of neutrons are due to the stronger values of the symmetry potentials when the 
momentum-dependence is turned off as shown in Fig. 4. 

C. Effects on isospin fractionation 

Isospin fractionation refers to the different isospin asymmetries in low and high density 
regions in nuclear reactions [34-36]. With symmetry energies increasing with density it is 
energetically more favorable to have the low density regions more neutron-rich while the high 
density regions more isospin symmetric. The strength of the isospin fractionation has been 
found to depend sensitively on the density-dependence of the symmetry energy [14,35,36]. 
As we have discussed in the previous subsection, the symmetry potential generally enhances 
the emission of neutrons while suppressing that of protons. A stiffer symmetry energy thus 
leads to a more neutron-rich low density gas of nucleons. Because of the charge and baryon 
number conservation, the remaining part (with local densities p > po/S) is consequently 
more isospin symmetric. It should be noted that there is a net conversion of neutrons to 
protons because of the more abundant production of Tr~ than tt"^ above the pion production 
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threshold. This effect, however, is very small even at a beam energy of 2 GeV/nucleon [7]. 
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FIG. 8. Isospin asymmetry of bound nucleons in central (upper window) and midcentral (lower 
window) collisions of ^^^Sn -)-i24 ^ beam energy of 400 MeV/nucleon. 

To first explore effects of the symmetry potential on isospin fractionation we study in Fig. 
8 the isospin asymmetry (neutron/proton ratio) of bound nucleons as a function of time for 
central (upper window) and midcentral (lower window) collisions of ^^'^Sn+^'^'^ Sn at a beam 
energy of 400 MeV/nucleon. First of all, from the lowest three curves (the three calculations 
use three different isoscalar potentials but the same symmetry potential of U^^^) in both 
windows, it is seen that the n/p ratio of bound nucleons (in the residue) is only slightly 
affected by the stiffness and the momentum-dependence of the isoscalar potential. They 
are calculated with the same symmetry potential U^'^^{p) but three different isoscalar 



19 



potentials, the MDYI, soft BKD and the stiff BKD. There is only a variation of about 1-3% 
in the (?^/p)p>po/8 ratio in the midcentral and central collisions. While it is very obvious 
that the (ri/p)p>po/8 ratio is much more sensitive to the variation of the symmetry potential. 
It is seen that the isospin asymmetry of bound nucleons increases by about 12% and 5% in 
central and midcentral collisions, respectively, going from the stiffer to the softer symmetry 
energy (from the U^^^ through MDI(O) to the MDI(l)). 

What are the effects of the momentum-dependence of the symmetry potential on the 
isospin fractionation? We study this question in Fig. 9 by comparing results obtained with 
and without the momentum-dependence of the symmetry potential. It is interesting to see 
that although the momentum-dependence of the symmetry potential increases the overall 
magnitude of isospin asymmetry of bound nucleons, relative effects of the symmetry energy 
remain about the same. In fact, in central collisions with the momentum- dependence of the 
symmetry potential effects of the symmetry energy become even slightly stronger. 
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FIG. 9. Effects of the momentum-dependence of the symmetry potential on isospin asym- 
metry of bound nucleons in central (upper window) and midcentral (lower window) collisions of 
^^^S'n +124 beam energy of 400 MeV/nucleon. 
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FIG. 10. Rapidity (upper window) and transverse momentum (lower window) distributions of 
isospin asymmetry of free nucleons with the three isoscalar potentials in the ^^^571+^^^ Sn reaction 
at a beam energy of 400 MeV/nucleon and an impact parameter of 5 fm. 



D. Combined constraints on density- and momentum-dependence of symmetry 

potentials 

Since the n/p ratio of bound or free nucleons are mainly affected by the symmetry 
potential unlike the inclusive observables which are strongly affected by isoscalar potentials, 
we concentrate in this subsection on the isospin asymmetry of free nucleons. Our main focus 
is to find observables that can be combined together to constrain simultaneously both the 
density and the momentum dependence of the symmetry potential. Specifically, we study 



22 



the isospin asymmetry of free nucleons as a function of rapidity and transverse momentum. 
First, to show that these observables are indeed not so sensitive to the isoscalar potential 
we show in Fig. 10 the isospin asymmetry of free nucleons as a function of rapidity (upper 
window) and transverse momentum (lower window) calculated using the same symmetry 
potential U^'^^{p) but the three different isoscalar potentials MDYI, soft BKD and stiff 
BKD. It is seen that there is a variation of less than 8% in the rapidity distribution of Sfree 
by changing from the MDYI to the soft BKD. The n/p ratio of free nucleons at midrapidity 
is shown as a function of transverse momentum in the lower window. For this observable the 
variation due to the change in isoscalar potential is even smaller at all transverse momenta. 
It should be mentioned that the overall rise of dN^/dNp at low transverse momenta is due 
to the Coulomb repulsion on protons. 

How can we constrain both the density- and momentum-dependence of the symmetry 
potential? Shown in the upper window of Fig. 11 are the isospin asymmetries of free nucle- 
ons calculated with and without the momentum-dependence of the symmetry potential as a 
function of rapidity. First of all, there is a shift downward when the momentum-dependence 
of the symmetry potential is included for both the soft (MDI(l)) and the stiffer (MDI(O)) 
symmetry energies. This is what we expected based on the comparisons of symmetry poten- 
tials in Fig. 4. Since the strength of the symmetry potential decreases with the increasing 
momentum, the reduction of Sfree is particularly larger at higher rapidities. It is also inter- 
esting to note that effects of the density-dependence of the symmetry potential, particularly 
stronger at midrapidity, arc about the same with or without the momentum-dependence 
of the symmetry potential. However, the value of Sf^ee with the soft symmetry energy but 
without the momentum-dependence of the symmetry potential (dotted line) is about the 
same as the one obtained using the momentum-dependent stiff symmetry energy MDI(O) 
(long dashed) at midrapidity, although they are quite different at forward and backward 
rapidities. 
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FIG. 11. Rapidity and transverse momentum distributions of isospin asymmetry of free nucle- 
ons with and without the momentum dependence of the symmetry potential in the ^^'^Sn +^^4 
reaction at a beam energy of 400 MeV/nucleon and an impact parameter of 5 fm. 

How do we separate them at midrapidities? One possible way is shown in the lower 
window of Fig. 11 where we compare the ratio of free neutrons to protons at midrapidity as a 
function of transverse momentum with the MDI(O) and the MDYl+U^y^^^^^p) interactions. 
It is seen that the two interactions are clearly separated at high transverse momenta. This is 
mainly because the symmetry potential U^^^^^^ (p) remains a constant while that embedded 
in the MDI(O) interaction decreases with increasing momenta as shown in Fig. 4. Thus 
by combining the measurements of the rapidity and transverse momentum distributions 
of the isospin asymmetry of free nucleons one can determine simultaneously the density- 
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and momentum-dependence of the symmetry potential. The density dependence of the 
symmetry energy can then be determined uniquely. We stress here that the separation at 
high transverse momenta is not due to the weaker isoscalar potential at high pt with the 
MDI(O) interaction as shown in Fig. 5. This is because the isoscalar potential acts equally on 
both neutrons and protons, the ratio dN^/dNp is thus insensitive to the EOS of symmetric 
nuclear matter as shown in Fig. 10. 




FIG. 12. Transverse momentum distributions of isospin asymmetry of free nucleons with and 
without the momentum dependence of the stiff (upper window) and soft (lower window) symmetry 
potential in the ^^^5n +^^4 reaction at a beam energy of 400 MeV/nucleon and an impact 
parameter of 5 fm. 

For a given density-dependent symmetry energy Egymip), how does the momentum- 
dependence of the symmetry potential affect the pt distribution of dNn/dNpl To answer this 
question we compare in Fig. 12 the dNn/dNp distributions with the stiff (upper window) 
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and soft (lower window) symmetry energy with (open and filled circles) and without (open 
and filled squares) the momentum-dependence of the symmetry potential. It is seen that 
without the momentum-dependence of the symmetry potential, the dNn/ dNp ratio turns to 
saturate above about pt > O.SGeV/c. With the momentum-dependence of the symmetry 
potential, however, the dNn/ dNp continue to decrease to high pt values. Again, this is a 
direct reflection of the decreasing symmetry potential with the increasing momentum. From 
both Fig. 11 and Fig. 12 one sees clearly the importance of measuring the isospin asymmetry 
of high Pt particles at midrapidity. 



E. Effects on rapidity and transverse momentum dependence of free nucleons in 

central collisions 

In the previous subsection, we concentrated on midcentral collisions at an impact pa- 
rameter of 5 fm. To examine the impact parameter dependence, we examine effects of the 
momentum-dependence of the symmetry potential in the ^^'^Sn -|-^^^ Sn reactions at 400 
MeV/nucleon and an impact parameter of 1 fm with the soft and stiff symmetry energy in 
Fig. 13 and Fig. 14, respectively. Comparing the results at an impact parameter of b=l 
fm in Fig. 13 and Fig. 14 with those at an impact parameter of b=5 fm in the previous 
subsection, it is seen that the main effects of the momentum-dependence of the symmetry 
potential do not change that much. This is mainly because for such a heavy reaction system 
as the ^^'^Sn -1-^^'^ Sn, the variation of impact parameter from 1 fm to 5 fm does not result in 
a very appreciable change in both the density and the isospin asymmetry in the participant 
region during the reaction. More quantitatively, with the soft symmetry energy MDI(l) the 
effects of the momentum-dependence of the symmetry potential at b=l fm are somewhat 
stronger than that at b=5 fm. However, with the stiff symmetry energy the impact param- 
eter dependence is less obvious. This is understandable from the strength of the symmetry 
potentials shown in Fig. 4. It is seen there that the slope of the symmetry potential depends 
more strongly on the density in the case of the soft symmetry energy (MDI(l)). Thus, as the 
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density increases due to the increasing centrality, the force due to the symmetry potential 
increases more strongly with the soft symmetry energy, leading to the larger effects shown 
in Fig. 13. 




FIG. 13. Rapidity (upper window) and transverse momentum (lower window) distributions 
of isospin asymmetry of free nucleons with and without the momentum dependence of the soft 
symmetry energy in the ^"^"^Sn +1^^ Sn reaction at a beam energy of 400 MeV/nucleon and an 
impact parameter of 1 fm. 
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FIG. 14. Rapidity (upper window) and transverse momentum (lower window) distributions 
of isospin asymmetry of free nucleons with and without the momentum dependence of the stiff 
symmetry energy in the ^^^Sn -j-^^*^ Sn reaction at a beam energy of 400 MeV/nucleon and an 
impact parameter of 1 fm. 



F. Effects on proton direct flow and neutron-proton differential flow 

How does the momentum-dependence of the symmetry potential affect the nuclear col- 
lective flow? First of all, it is important to know what collective observables are sensitive 
to the symmetry potential with or without the momentum- dependence. Since symmetry 
potential is small compared to the isoscalar potential, moreover, at high energies baryon- 
baryon collisions dominate the reaction dynamics, one has to use delicate observables to 
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extract information about the symmetry potential. To illustrate this point, we show in Fig. 
15 the standard transverse flow analysis for protons in the reaction of ^^"^Sn +^^^ 5*^ at a 
beam energy of 400 MeV/nucleon and an impact parameter of 5 fm. It is seen that with 
(upper window) or without (lower window) the momentum-dependence of the symmetry 
potential for both the soft and stiff symmetry energies, the proton transverse flow is about 
the same. 




FIG. 15. Proton transverse flow analysis with (upper window) and without (lower window) 
the momentum-dependence of the symmetry potential in the reaction of ^^"^Sn +1^4 ^ beam 

energy of 400 MeV/nucleon and an impact parameter of 5 fm. 
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FIG. 16. Neutron-proton differential flow analysis with (upper window) and without (lower 
window) the momentum-dependence of the symmetry potential in the reaction of ^^'^Sn +124 
at a beam energy of 400 MeV/nucleon and an impact parameter of 5 fm. 

We now turn to the neutron-proton differential flow 

N{y) 

F:_^{y)^Y.iP>^)/N{y), (16) 

■t=i 

where Wi = 1(— 1) for neutrons (protons) and N{y) is the total number of free nucleons at 
rapidity y. The differential flow combines constructively effects of the symmetry potential on 
the isospin fractionation and the collective flow. It has the advantage of maximizing effects of 
the symmetry potential while minimizing effects of the isoscalar potential [14] . Shown in Fig. 
16 is the neutron-proton differential flow with (upper window) and without (lower window) 
the momentum-dependence of the symmetry potential. First, it is seen that in both cases, 
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the neutron-proton differential flow is rather sensitive to the symmetry energy compared to 
the standard proton transverse flow analysis shown in Fig. 15. Secondly, it is seen that 
the neutron-proton differential flow is stronger without the momentum-dependence of the 
symmetry potential. This is because the momentum-independent symmetry potential makes 
not only more neutrons free but also gives them higher transverse momenta in the reaction 
plane because of its higher magnitude and steeper slopes compared to the momentum- 
dependent symmetry potential embedded in eq. 2. 

Is the neutron- proton differential flow still sensitive to the variation of Egymip) when the 
momentum-dependent symmetry potential is used? The answer is a definitely yes as shown 
in the upper window of Fig. 16. In fact, the slope of the differential flow at y=0 

changes from about -1.9 to 21.4 MeV/c using eq. 2 by changing from the MDI(l) to the 
MDl(O) parameter set, while it changes from about 26 to 39 MeV/c with the momentum- 
independent symmetry potentials. The net change due to the variation of the Egymip) is 
thus actually larger with the momentum-dependent symmetry potentials. Thus the neutron- 
proton differential flow remains a sensitive probe to the density dependence of the symmetry 
energy. 

IV. SUMMARY 

In summary, using an isospin- and momentum-dependent transport model we studied ef- 
fects of the momentum-dependent symmetry potential on several experimental observables 
in heavy-ion collisions induced by neutron-rich nuclei at RIA energies. The momentum- 
dependence of the symmetry potential is found to play an important role in the isospin 
fractionation, particle emission and neutron-proton differential flow. In particular, high 
Pt neutron to proton ratio at midrapidity is very sensitive to the momentum-dependence 
of the symmetry potential. However, the momentum- dependence of the symmetry poten- 
tial has httle effects on the overall reaction dynamics and the directed nuclear flow. To 
uniquely determine the Esym{p) one needs to determine simultaneously both the density- 
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and momentum-dependence of the symmetry potential. For this purpose, it is necessary to 
combine several experimental observables, such as, the rapidity and transverse momentum 
dependence of the free neutron to proton ratio, and the neutron-proton differential flow. 
Comparing experimental data with predictions based on transport models including the 
momentum-dependent symmetry potential is thus crucial for investigating accurately the 
EOS of dense neutron-rich matter. 
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